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Abstract. The combination of a thin disk and a narrowly-collimated jet is a typical 
structure that is observed in the vicinity of a massive object such as AGN, black hole or 
YSO. Despite a large variety of their scales and possible diversity of involved processes, 
a simple and universal principle dictates the geometric similarity of the structure; we 
show that the singularity at the origin (r = 0) of the Keplerian rotation (Ve cx r^^/^) 
is the determinant of the structure. The coUimation of jet is the consequence of the 
alignment — so-called Beltrami condition — of the flow velocity and the "generalized 
vorticity" that appears as an axle penetrating the disk (the vorticity is generalized 
to combine with magnetic field as well as to subtract the friction force causing the 
accretion) . Typical distributions of the density and flow velocity are delineated by a 
similarity solution of the simplified version of the model. 
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1. Introduction 

An accretion disk often combines with spindle-like jet of ejecting gas, and constitutes a 
typical structure that accompanies a massive object of various scales, ranging from young 
stars to galactic nuclei |20] . The mechanism that rules each part of different systems 
might be not universal. Since early 70s after the discovery of radio galaxies and quasars 
(see e.g. [3] and references therein) the main evidence from detecting jets in different 
classes of astrophysical systems which are observed to produce coUimated jets near the 
massive central object prove the direct association with an accretion disk, although 
possibly reflecting different accretion regimes; the opposite is not true in some objects 
for which the accretion disks do not require coUimated jets (viscous transport /disk winds 
may play the similar role in the energy balance) [TJ |28l [13]. The macroscopic disk-jet 
geometry, however, bears a marked similarity despite the huge variety of the scaling 
parameters such as Lorentz factor, Reynolds number, Lundquist number, ionization 
fractions, etc. [TJ El EH [291 [131 HE]. The search for a general principle that dictates such 
similarity provided the stimulus for this paper. 

Before formulating a model of the global structure, we start by reviewing the basic 
properties of disk-jet systems. In the disk region, the transport processes of mass, 
momentum and energy depend strongly on the scaling parameters; many different 
mechanisms have been proposed and examined carefully. The classical (coUisional) 
processes are evidently insufficient to account for the accretion rate, thus turbulent 
transports (involving magnetic perturbations) must be invoked |2]. Winds may also 
remove the angular momentum from the disk [211 [38l [TO] . The connection of the disk 
and the jet is more complicated. While it is evident that the mass and energy of the 
jet are fed by the accreting flow, the mechanism and process of mass/energy transfer 
are still not clear. It is believed that the major constituent of jets is the material of 
an accretion disk surrounding the central object. Although for the fastest outflows the 
contributions to the total mass flux may come from outer regions as well |36]. In AGN 
one may think of taking some energy from the central black hole. Livio [28] summarizes 
the conjectures as follows: (i) powerful jets are produced by systems in which on top 
of an accretion disk threaded by a vertical field, there exists an additional source of 
energy/wind, possibly associated with the central object (for example, stellar wind from 
porotostar may accelerate YSO jets, as estimated by Refs [III [15l [IS]); (ii) launching of 
an outfiow from an accretion disk requires a hot corona or a supply from some additional 
source of energy (see also [HI [261 [36]); (iii) extensive hot atmosphere around the compact 
object can provide additional acceleration. Magnetic fields are considered to play an 
important role in defining the local accretion [3 [HI [HI [221 ESI SS] • When magnetic field 
is advected inwards by accreting material or/and generated locally by some mechanism, 
i.e. if there is a bundle of magnetic field lines which is narrowest near the origin (central 
object) and broader afar off, the centrifugal force due to rotation may boost the jet along 
the magnetic field lines up to a super-Alfvenic speed [H [H [5l [12l [25l [271 [1] • In the case of 
AGN, there is an alternative idea suggested by Blandford & Znajek [7] based on electro- 
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dynamical processes extracting energy from a rotating black hole. Extra-galactic radio 
jets might be accelerated by highly disorganized magnetic fields that are strong enough 
to dominate the dynamics until the terminal Lorentz factor is reached [I9] . Following 
the twin-exhaust model by Blandford & Rees [6] , the coUimation under this scenario is 
provided by the stratified thermal pressure from an external medium. The acceleration 
efficiency then depends on the pressure gradient of the medium. In addition to the 
energetics to account for the acceleration of ejecting fiow, we have to explain how the 
streamlines/magnetic-field lines change the topology through the disk-jet connection jlQ] 
(see also [391 E]). 

Despite the diversity and complexity of holistic processes, there must yet be a simple 
and universal principle that determines the geometric similarity of disk-jet compositions. 
In the present paper we will show that the collimated structure of jet is a natural 
consequence of the alignment of the velocity and the "generalized vorticity" — here the 
notion of vorticity will be generalized to combine with magnetic field (or, electromagnetic 
vorticity) as well as to include the effective friction causing the accretion. In a Keplerian 
thin disk, the vorticity becomes a vertical vector with the norm oc r~^/^ (r is the radius 
from the center of the disk), which appears as a spindle of a thin disk. Then, alignment 
is the only recourse in avoiding singularity of force near the axis of Keplerian rotation. 
We are not going to discuss the "process" that produces the global structure; we will, 
however, elucidate a "necessary condition" imposed on the structure that can live long. 

To demonstrate how the alignment condition arises and how it determines the 
singular structure of a thin disk and narrowly- collimated jet, we invoke the simplest 
(minimum) model of magnetohydrodynamics (Sec. [2]). The alignment condition (so- 
called Beltrami relation) will be derived in Sec. El In Sec. HI we will formulate a system 
of equations which gives general Beltrami structures satisfying the alignment condition. 
In Sec.[5l we will give an analytic solution to a simplified version of the equations by a 
similarity-solution method. 

2. Momentum Equation 

We start by formulating a magnetohydrodynamic (MHD) model. Let P = pV denote 
the momentum density, where p is the mass density and V is the (ion) fiow velocity. 
The momentum balance equations for the ion and electron fiuids are 



where B is the magnetic field, (p and A are the the electromagnetic potentials, (j) is the 
gravity potential, pi and pe are the ion and electron pressures, and z/ is the (effective) 
friction coefficient. We are assuming singly charged ions, for simplicity. The variables 
are normalized as follows: We choose a representative fiow velocity Vq and a mass density 



dtP + V ■ {pV ® V) 

= e-^p{V xB- 
- pV(0 + ifi) 



dtA) 

- Vp, - uP, 



(1) 
(2) 
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Po in the disk, and normalize V and p by these units. The energy densities |Sp/87r , 
cif, p(j), Pi and Pe are normahzed by the unit kinetic energy density 

= (3) 

The independent variables (coordinate x and time t) are normalized by the system size 
Lq and the corresponding transit time To = Lq/Vq. The scale parameter e is defined by 

s. 

where Si = mc/ y/Ane^pQ is the ion inertia length (m: ion mass, e: elementary charge). 

In equation ([2]), we are neglecting the inertia and the gravitation of electrons, as 
well as the electric resistivity. The electron velocity is given by 

Ve = V - ep-^V X B. 

The system ([I])-© includes a finite dissipation by the effective friction force on the (ion) 
momentum, not by conventional viscosities or resistivity. This is for the simplicity of 
analysis. We do need a finite dissipation to allow accretion, but the detail mechanism 
of dissipation is not essential in the scope of our arguments. 

In this paper, we consider stationary solutions, so we set dt = 0. To study the 
macroscopic structure of an astronomical system, we may assume e ^ 1. Let us expand 
B as 

B = + eB(i) + ■ ■ ■ , 

where the magnitude of each S^") is of order unity. All other independent and dependent 
variables are assumed to be at most of order unity. Then, equation ([2]) reads as 

= e-^pV X 

+ p[Vx - p~\V X B(°)) X 

- pVip + Vpe + 0{e). 

The term of order demands 

= pP, (4) 

where /i is a certain scaler function, which is the reciprocal Alfven Much number. 
Operating divergence on both sides of equation (j4]), we find V ■ (pP) = P ■ V p = 0. 
From order e*^ terms, we obtain 

V X = p~\V X B(°)) X + V<^ - p-^Vpe. 
Using these relations in equation ([1]), we obtain 

V ■ (py ® V) = [V X (pP)] X (pP) - pWcj) -Vp- uP, (5) 
where p = pi + pe- 

In order to derive a term that balances with the friction term, we decompose the 
"inertia term" [the left-hand side of equation ([5])] as follows: we first write 

p = PlP2 (6) 
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(p2 will be determined as a function of u in equation (|TT]) ). and denote 

Pi = p,V, P2 = P2V. (7) 

Using these variables, we may write 

V -{pV ®V) = V -{Pi® P2) 

= (V • Pl)P2 + (Pi ■ V)P2. (8) 

In the conventional formulation of fluid mechanics, we choose P2 = 1 and pi = p. Then, 

V ■ (pV ® V) = (V ■ P)V + p{V ■ V)V. Combining with dtP = pdtV + V dtp, and 
using the mass conservation law dtp + V ■ P = 0, we observe that the left-hand side of 
equation ([5]) reduces into the standard inertia term p(V ■ V)V (here we are considering 
the steady state with dtV = 0). In the preset analysis, however, we choose a different 
separation of the inertia term to match the term (V ■ Pi)P2 with the friction term —uP; 
this strategy will be explained in the next section. 

Multiplying (P2/P1) on both sides of equation ([5]) , we obtain 

P2xn = ^vp^ + plv{<p + h) 

+ P2Z/P2 + ^(V-Pl)P2, (9) 

Pi 

where 

= V X P2 - PP2V X (p p) (10) 

= V X P2 - PP2V X 
is a generalized vorticity, and h is the enthalpy density (Vh = p^^Vp). 

3. Disk-Jet Structure 

Here we assume toroidal symmetry {dg = in the r-6-z coordinates). We consider a 
massive central object (a singularity at the origin), which yields (p = —MG/r (we may 
neglect the mass in the disk in evaluating 0). Then, the flow velocity is V ~ Vgeo with 
the Keplerian velocity Vg oc r~^/^ in the disk region. The vorticity is V x V = Vtze^ 
with VLz oc r~^/^. The momentum is strongly localized in the thin disk, and the vorticity 
diverges near the axis (in this macroscopic view). This singular configuration allows 
only a special geometric structure to emerge. By equation (|9]), following conclusions are 
readily deducible: 

(i) In the disk, a radial flow (which is much smaller than Ve) is caused by the friction. 
The friction force P2V P2 is primarily in the azimuthal (toroidal) direction, which may 
be balanced with the term (P2/P1) (V • Pi)P2 that has been extracted from the inertia 
term V ■ {pV ® V); see equation ([8]). Using the steady-state mass conservation law 

V ■ P = 0, we observe 

pr V ■ Pi = pr'v ■ (P/P2) 

= P2V- Vp2^ = - V- Vl0gp2. 
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Hence, the balance of the friction and the partial inertia term demands 

V-V\ogp2 = iy, (11) 

which determines the parameter p2. The remaining part pi of the density is, then, 
determined by the mass conservation law: From 

V ■ P = V ■ (pV) 

= P2V ■ Vpi + piV ■ Vp2 + P1P2V ■ V 

= 

and equation (fTTj) . we obtain a relation 

V- Vlogpi = -V- V-z/. (12) 

(ii) After balancing the third and fourth terms in equation ([9]) , the remaining 
terms must not have an azimuthal (toroidal) component. In fact, the right-hand side 
(gradients) have only poloidal {r-z plane) components, and hence, the left-hand side 
cannot have a toroidal component. 

(iii) In the vicinity of the axis threading the central object, the flow V{= P2/P2) 
must align to the generalized vorticity Q (that is dominated by V x V oc r^'^/^e^), 
producing the coUimated jet structure, to minimize the Coriolis force V xfl: Otherwise, 
the remaining potential forces (gradients of potentials which have only poloidal 
components) cannot balance with the rotational (toroidal) Coriolis force. The alignment 
condition, which we call Beltrami condition, reads as 

n = AP, (13) 

where A is a certain scalar function. The pure electromagnetic Beltrami condition is the 
well-known force-free condition V x B = \B [HI 1301 112] • Equation f|T3|) demands the 
alignment of the generalized vorticity and momentum [32j . 

(iv) When the Beltrami condition eliminates the left-hand side of equation ([9]), the 
remaining potential forces must balance and achieve the Bernoulli condition [32] which 
reads as 

^VP| + V((/) + /i) = V {U'^ + (t) + h^+ l^'Vlogp2 

= 0. (14) 

The system of determining equations is summarized as follows: By equation (iTTl) . we 
determine the "artificial ingredient" p2 for a given v. This equation involves V = P/p 
which is governed by the Beltrami equation f lT3|) . After determining V and p2, we 
can solve the Bernoulli equation (HM to determine the enthalpy h (the gravitational 
potential is approximated by = —MG/r). 

4. Beltrami Vortex Structure 

4-1. General two-dimensional structure 

In this section, we rewrite the Beltrami-Bernoulli equations (fTT]) -( [T^ in a more 
manageable form by invoking the Clebsch parameterization [45]. In an axisymmetric 
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geometry, the divergence-free vector P may be parameterized as 

where I = prVg. Both ip and / do not depend on 9. Since P ■ Vip 
(contours) of if) are the streamhnes of P (or those oiV = P / p). 
Using the expression (fTSjl in equation (ITT!) yields 



(15) 

0, the level sets 



1 
P 



P- Vlogps = -Vlogpa X VV^- 



P 
1 

rp 



{logp2,V}, 



(16) 



where {«.,&} = {drb){dza) — {dra){dzh). For a given set of P, p and z/, we can solve 
equation (|T6l) to determine p2, as well as pi = p/p2 which is consistent to equation (fT2l) . 

Let us rewrite the momentum equation ([9]) using the Clebsch parameterization f|T5|) . 
We observe 



17 = V X 



+ 



- mV X (/iP) 



P>2U^ + VV^- Vpr'-;Up2Vp 



v(pr'/) -PP2V(p/)] X v^^, 



(17) 



where 



Cip = rdr{r ^drip) + d'^ijj. 
As mentioned above, P2 x may not have a toroidal component, i.e.. 



V (pr'/) - PP2V (p/)] X = 0, 

which is equivalent to the existence of a scalar function A such that 
AV^ = V (pi'l) - PP2V (pJ) 

= (pr' - p'p2)v/ + jvpr' - pp2/vp. 

The poloidal component of the momentum equation (Q reads 
{pY' - PV2) (C^) 

■ (vp^^ - pp2Vp)] Vip + IV I 

= -yV [r-Vr^(|V^p + /^); 
-rV^V(0 + /i). 



(18) 



(19) 



We have to determine a self-consistent set of functions J, /i, pi, p2, A and p (the 
friction coefficient v and the gravitational potential = —MG/r are given functions). 
Our system of equations are (fT6ll . ( ITSll and ( !T9ll . which are simultaneous nonlinear partial 
differential equations including hyperbolic parts. Since P ■ Vp = [see equation (jl])], 
we may assume that p = p{ip), and give this function as a "Cauchy data" (an arbitrary 
function of ip that is constant along the characteristics). Another function (for example 
Pi) may be arbitrarily prescribed to leave five unknown functions to be determined by 
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five equations of tlie system. However, the fiyperbolic parts must be integrated by some 
Caucliy data, not by boundary values. We note tliat tfie factor (pj"^ — /^^P2) multiplying 
the elliptic operator £ in equation ( |T9l) may cause the Alfven singularity. 

To proceed with analytic calculations, we will consider simplified systems in which 
we can integrate the hyperbohc part of the equations easily. 



4-.2. Unmagnetized Beltrami flow 

Here, we consider unmagnetized (or super-Alfvenic; /i ~ 0) Beltrami-Bernoulli flows. 
With /i = 0, equation ( ITSi) reduces to 

AW = V (pr'z) , (20) 

which implies pi^I = hi'ip) (Cauchy data) and A = \{tp) = l2{'4') (we denote 
fi^p) = df{ip)ldip). We notice that equation (l20l) is just the toroidal (azimuthal) 
component of the Beltrami condition ( fT3l) . Applying the Beltrami condition also to 
the poloidal component of the momentum equation, we can simplify equation (fT9|) : with 
/i = and / = pil2{ip), we obtain 

Cij - ■ Vlogpi = -p?/^(^)/2(^). (21) 

The Beltrami condition has decoupled the gradient forces [the right-hand side of 
equation (1191) ] from the momentum equation, which must balance separately: This is 
the Bernoulli condition ([TH) which now reads as 



Vh= - V 
1 



+ /2)Vlogp2. (22) 



Now the determining equations are much simplified — for a given distribution pi 
and Cauchy data hii^), we may solve equation (!?!]) to determine ip- Then, for a given 
z/, equation (|T6|) yields p2. Finally, the enthalpy h is determined by equation fl22|) . 



4-3. A simple magnetized Beltrami flow 
If we assume 

Pi = PiW, (23) 

(i.e., V ■ Vpi = and thus V ■ P2 = 0), we may rather easily include a magnetic field 
into the solution (p = fi{ip) 7^ 0). By equation (IT^ . this assumption implies that the 
compressibility is only by the frictional deceleration (then, = V ■ P = V ■ (P1P2) = 
V-Ps). 

Under the assumption (1231) . the second and third terms of equation ( !T8|) parallel 
Vijj. Therefore, V/ must align to Vip, implying / = In this case, A is given by 

A = (pr' - MV2)/' - (prV'i + P2/ip )/• (24) 
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Figure 1. The momentum field (contours of that describe the streamlines of the 
poloidal component of P) of the similarity solution (with D = 1, p = 1, J = 0.1 and 

9 = 1)- 



On the right-hand side of equation fl24p . only p2 is not a function of ip. The Behrami 
condition demands vanishing of the left-hand side of equation f|T9|) . which, with a Cauchy 
data I{ip) and the A of equation ( l24|) . reads as 

(pr' - /iV2)/:^ - {pfp'i + P2/i/i')|v^r = (25) 

This equation (governing ip) is coupled with equation ( fT6ll through p2. The Bernoulli 
condition (1221) can be separately solved to determine h. 



5. Analytic Similarity Solution 

5.1. A similarity solution modeling disk- jet structure 

In this section, we construct a similarity solution of the un-magnetized (p = 0) model 
( |2T|) . which describes a fundamental disk-jet structure. We define 

r = - (r > 0), (26) 
r 

and an orthogonal variable (Vr ■ Vcr = 0) 



0" 



(27) 



We consider such that 

= ^(r) = -JtP - Dt-'^, (28) 

where J and p (D and q) are positive constants, which control the strength of the jet 
(disk) flow. As shown in Fig.[Tl this ip has a disk-jet-like geometry. 

The level sets of ip (hence, those of r) are the streamlines of P. On the other hand, 
a serves as the coordinate directed parallel to the streamlines. We assume that pi is 
written as 

Pi(r,a)=px(r)p||(a), (29) 
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Figure 2. The distribution of p± of the similarity solution (with D 
J = 0.1 and q = 1). Contour curves are given in the log scale of p±. 



1, P = 1, 



and, then, logpi = log p^(r) + logp||(cr). 

Let us see how the stream function of equation (128|) satisfies equations (12T!) , (l22l) 
and (fT6ll . i.e., we determine all other fields hi'^P), pi, P2, ^ and h that allow this to 
be the solution. For arbitrary fir) and gir)^ we observe 

-C/ = ^ [(^' + 1)/" + 3r/' 



Hence, the left-had side of equation (pT]) is (denoting gij) = logp_L(r)) 
Cip — Vip ■ V log pi 

= \ \{t^ + + Sri;' - (r^ + l)g'i;' 



(30) 



For this quantity to balance with the right-hand side of equation f l2T|) . which has no 
explicit dependence on r~^, both sides must be zero, i.e., we have to set 12(4') = (the 
implication of this simple condition will be discussed later). Then, equation fl^T]) reduces 
to 

(r^ + 1)^" + 3t^' - (r^ + l)g'^' = 0. (31) 
We note that the Beltrami condition ( 13T1) is freed from P\\{(t). This fact merits in solving 



equation (11611 : see equation ([3^ 

For the specific form (|28l) of tp, we have to determine an appropriate g = log p± to 
satisfy equation fl3Tl) . i.e.. 



_ J p {p - 1)tP+'' + D q {q + 1] 
~ J p t P+i+'^ - D q T 



+ 



3r 



r2 + 1 



(32) 
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Figure 3. The density p in tlie similarity solution (with D ~ 1, p = 1, J = 0.1 and 
q — 1). We assume p\\{<7) ~ z^^^^ and p2 ~ A levelset surface of p is shown in 

the doiman r < 5 and \z\ < 5. 



Solving equation (l32l). we obtain 

^ = logp± = log + - log(r + 1), 

and, thus, 

= • (33) 

In Fig. [21 we show the profile of p±(r). 

5.2. Bernoulli relation in the disk region 

As mentioned above, this solution assumes /2(V') (= A) = 0, and hence, I2 = p2fVe must 
uniformly distribute. In the disk region (the vicinity of 2; = 0), we may approximate 
Ve ~ \J MG I r (Keplerian velocity). Hence, p2 oc r~^/^. In Fig.|3l we show the profile of 
P = P1P2 for the case of pi oc pj_ (i.e., p\\ =constant). 
For p = P1P2 = p\\p±r~^^'^, equation f|T6l) reads as 
—dib/dz o,, J.-5/2 
2p||P± p\\{r) 

along each streamline in the disk region. For a given i/, we can solve equation (1341) for 
P|j to determine the density profile. 

In the disk region the Bernoulli relation ( IT4l) accounts as follows: by V • P = 0, we 
have Pr = pVr oc r~^ . If P||(cr) = constant, for example, p oc r~^/^ (evaluated along a 
streamline in the disk region). Then, we have = Ko'""^^^ with a (negative) constant 
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Vro- Combining the azimuthal velocity Vg = Vqqt (which must be shghtly smaller 
than the Keplerian velocity J MG/r)^ we obtain 



By p2 oc r~^/^, we obtain dr{logp2) = — (l/2)r^^. Hence, the Bernoulli relation ([1 
demands 

d,h = {V^ - MG)r-\ 

which yields h = {MG — V"g^)r~^. In this estimate, all components of the energy density 
(gravitational potential 0, kinetic energy and enthalpy h) have a similar profile 

(oc r~^). 

5. 3. Bernoulli relation in the jet region 

In the jet region (vicinity of r = 0), the streamlines (contours of r = 2;/r) are almost 
vertical, and we may approximate a ^ z. 

Let us first estimate p2 using equation (ITBl) . which is approximated, in the jet region, 

by 

rpu = {\ogp2,^p} 

fti (drip) {dAogp2) 

= JprP+'^d,i\ogp2), (35) 

which shows that p2 is an increasing function of \z\. Using p = P\\{o')p±{t)p2, we 
integrate equation fl35|) along the streamline (r =constant, a ~ 2;): 

P2 \P2j JP 

With this p2{z), we may estimate the toroidal (azimuthal) component of the velocity: 
Ve = l2/{rp2) = {I2T)/ {ZP2), where I2 and r are constant (the latter is constant along 
each streamline). We find that the kinetic energy Vg/2 of the azimuthal velocity 
decreases as a function of \z\ (both by the geometric expansion factor z'"^ and the 
friction damping effect P2^^). The steep gradient of the corresponding hydrodynamic 
pressure yields a strong boost near the foot point (z ~ 0). 

The poloidal component of the kinetic energy is estimated as follows: We may 
approximate 

' Jp- 



2 ^ I - I = — r^^^ Pdr)p\\{z)z^dz. (36) 



2p2^2 \^ ^ j.p+1 



(37) 



2p2;24 

Here, the vertical distribution of the density p = p±{j)p\\{a)p2 is primarily dominated 
by pii(cr) ^ p\\{z). 



Beltrami structures in disk-jet system 



13 



At long distance from the origin, the jet has a natural similarity property. For 
simplicity, let us ignore the effect of the friction {u = 0), and assume p2 = 1. Then, 
Pll oc yields {V^ + V^)/2 oc ^ which may balance with the gravitational 

potential energy = —MG\z\~^ . Note, that the azimuthal component of the kinetic 
energy disappears at large scale [Vq oc The Bernoulli condition ( IT4l) gives h that 

also has a similar distribution of oc \z\^^ . 

In Fig. [3], we show the profile of p = p±P||P2 with py oc (jet region) and 

P2 oc r^^/^ (disk region). 

6. Summary and Concluding Remarks 

We have shown that the combination of a thin disk and narrowly- collimated jet is 
the unique structure that is amenable to the singularity of the Keplerian vorticity; 
the Beltrami condition, imposing the alignment of flow and generalized vorticity, 
characterizes such a geometry. Here the conventional vorticity is generalized to combine 
with magnetic field (or, electromagnetic vorticity) as well as to subtract the friction 
force causing the accretion. 

In SectionEl we have given an analytic solution in which the generalized vorticity 
is purely kinematic (fi = V x P2 with the momentum P2 that is modified by the 
friction effect). The principal force that ejects the jet is, then, the hydrodynamic 
pressure dominated by Vq /2. Additional magnetic force may contribute to jet 
acceleration if the self-consistently generated large-scale magnetic field is sufficiently 
large [iHl EHl [371 HD 132] • Such structures can be described by the generalized model of 
Subsection l4.1[ 

The similarity solution has a singularity at the origin (where the model gravity 
(j) = —MGr~^ is singular), which disconnects the disk part (parameterized by D and 
q) and the jet part (parameterized by J and p). To "connect" both subsystems (or, 
to relate D, q, J, and p), we need a singular perturbation that dominates the small- 
scale hierarchy |44j . The connection point must switch the topology of the flow: In the 
poloidal cross section, the streamlines of the mass flow, which connect the accreting 
inflow and jet's outflow, describe hyperbolic curves converging to the radial and the 
vertical axes. On the other hand, the vorticity lines (or generalized magnetic field 
lines) thread the disk. This topological difference of these two vector fields demands 
"decoupling" of them. Instead of disconnecting them by a singularity, we will have to 
consider a small-scale structure in which the topological switch can occur [40j . The Hall 
effect (scaled by e) or the viscosity/resistivity (scaled by reciprocal Reynolds number) 
yields a singular perturbation (then, the generalized magneto- Bernoulli mechanism 
[331 [M] may effectively accelerate the jet-flow). In a weakly ionized plasma, the 
electron equation ([2]) is modified so that the Hall effect is magnified by the ratio of 
the neutral and electron densities; the ambipoler diffusion effect also yields a higher- 
order perturbation [23]. The mechanism of singular perturbation and the local structure 
of the disk-jet connection point may differ depending on the plasma condition near the 
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central object: in AGN, the plasma is fully ionized but rather collisional (we also need 
a relativistic equation of state with possible presence of pairs), and, in YSO, the plasma 
is partially ionized; the relevant dissipation mechanisms determine the scale hierarchy 
in the vicinity of the singularity. The study of above concrete cases was not the scope of 
present paper and will be considered elsewhere. Our goal was to show how the alignment 
of flow and generalized vorticity condition arises and how it determines the singular 
structure of a thin disk and narrow ly-collimated jet invoking the simplest (minimum) 
model of magnetohydrodynamics. 
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